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ABSTRACT 



NASA's Small Explorer Mission GEMS (Gravity and Extreme Magnetism 
SMEX), scheduled for launch in 2014, will have the sensitivity to detect and 
measure the linear polarization properties of the 0.5 keV and 2-10 keV X-ray 
emission of a considerable number of galactic and extragalactic sources. The 
prospect of sensitive X-ray polarimetry justifies a closer look at the polarization 
properties of the basic emission mechanisms. In this paper, we present analyt- 
ical and numerical calculations of the linear polarization properties of inverse 
Compton scattered radiation. We describe a generally applicable formalism that 
can be used to numerically compute the polarization properties in the Thomson 
and Klein-Nishina regimes. We use the code to perform for the first time a de- 
tailed comparison of numerical results and the earlier analytical results derived 
by Bonometto et al. (1970) for scatterings in the Thomson regime. Furthermore, 
we use the numerical formalism to scrutinize the polarization properties of syn- 
chrotron self-Compton emission, and of inverse Compton radiation emitted in 
the Klein-Nishina regime. We conclude with a discussion of the scientific poten- 
tial of future GEMS observations of blazars. The GEMS mission will be able 
to confirm the synchrotron origin of the low-energy emission component from 
high frequency peaked BL Lac objects. Furthermore, the observations have the 
potential to decide between a synchrotron self-Compton and external-Compton 
origin of the high-energy emission component from flat spectrum radio quasars 
and low frequency peaked BL Lac objects. 

Subject headings: Polarization, Radiation mechanisms: non-thermal, Relativistic 
processes, Scattering, X-rays: galaxies, Gamma rays: galaxies, Galaxies: active, 
Galaxies: jets, BL Lacertae objects: general 
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Introduction 



Although th e scientific potential of X-ray po 



Rees 



1975 



Lightman fc Shapiro 



arim e try observations has 



1975 



1976 



Meszaros et al 



ong 



been 



1988|), only 



appreciated (e.g. 

a single dedicate d X-ray polarimete r mission has been flown so far: the OSO-8 mission 



launched in 1975 ( jNovick et al. 



1977| ) . The mission was able to detect polarization degrees of 
a few percent for X-ray bright sources with Crab-like fluxes. This sensitivity was sufficient to 
detect the polarization of one cosmic source. The observations of the Crab nebula revealed 



evide nce for linear polarization at 2.6 keV and 5.2 keV at a level of ~20% ( jWeisskopf et al. 



19781 ). More recently, the I NTEGRAL satel 



Crab Nebula and Cyg X-l ( IDean et al. 



2008 



i te was used to study the polarization of the 



Forot et al. 



2008 



Laurent et al 



201 lh . Even 



though the observations suggest polarization degrees of >40% in the hard X-ray/7-ray 
regime, the systematic uncertainties on these results are large. NASA plans to launch 
the Small Explorer Mission GEMS (Gravity and Extreme Magnetism SMEX) in 2014, a 
dedicated mission for X-ray polarimetry in the 2-10 keV energy range, with the sensitivity 
to detec t polarization deg rees down to ~ 1% even for weak sources with a flux of a few 



mCrab ( iBlack et al. 



2010). A student e xperiment wi 



capabilities into sub-keV energy range ( IKaaret et al. 



exte nd the polarimetric measurement 



20091 ). GEMS will be able to measure 



the polarization properties of dozens of galactic and extragalactic sources and to establish 
X-ray polarimetry as an observational discipline. The main science drivers of GEMS are 
the measurement of the orientation and inclination of the inner accretion disks and the 
masses and spins of Binary Black Hole systems (BBHs), the detection of plasma physics 
effects in neutron star magnetospheres, the study of the locale and properties of the particle 



jets c 


f Active Galac 


ic Nuclei (AGNs) 


see 


Lei. Dean. Hills et al. 


1997; 


WeisskoDf et al. 


2006; 


Bellazini et al. 


2010; 


Krawczvnski 


201 


L). 



-4- 



The upcoming launch of the GEMS mission motivates a close look at the polarization 
properties of the most common X-ray emission processes. Whereas the polarization 



emission are well studied (e.g . 



Biornsson Sz Blumentha] 



Ruszkowski fc Begelman 



1982 



2002 



Ginzbure 


& Svrovatskii 


1965 




Dolan 


1967 


Biornsson 


1985: 


Rvbicki & Lightman 


1986 



Laing 


198( 


]; 


Begelman 


1993 



and references therein), a similar in-depth understanding of 
the polarization properties of inverse Compton emission is still lacking. In the following, 
"Inverse Compton" denotes the emission from a single up-scattering of low-frequency target 
photons by relativistic electrons. Prime examples of inverse Compton emission include the 
high-energy component of AGNs and possibly 7-ray emission from Gamma-Ray Bursts. 
This paper focuses on the polarization of inverse Compton emission and the application 
of the results to AGN jets. We will not discuss the polarization of "Comptonized" 
emission (e. g. from AGN coronae), which underwent many Compton and inverse Compton 
scatterings before leaving the source. Inverse Compton scatterings tend to reduce the 
polarization degree of emission; multiple Compton scatterings do so even more. 

The polarization of inverse Compton emission usually arises from averaging the results 
from elemental scattering processes over a range of electron energies and propagation 
directions. Highly relativistic electrons with Lorentz factors 7 ^> 1 scatter photons into a 
narrow cone centered on their propagation direction with an opening angle of 2/7. The 
inverse Compton emission from a relativistic electron plasma thus comes from electrons 
with velocities aligned to within ~ 7 _1 with the line of sight. In this paper, we discuss the 
case that the angular distribution of the electrons contributing to the observed emission 
does not deviate appreciably from isotropy in the rest frame of the emitting plasma, called 
the plasma frame (PF) in the following. For large 7-values most electron populations will 
satisfy this criterium. 
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The paper of iBonometto. Cazzola. &: Saggionl (119701 ) (called BCS in the following) 
reports on the polarization of the inverse Compton emission on the basis of a quantum 
mechanical scattering calculation. The study makes a series of approximations and is 
limited to highly relativistic electrons (7 ^> 1) and scatterings in the Thomson regime 
•y(hu/m e c 2 ) <C 1 where oj denotes the frequency of the target photons. The paper presents 
equations to compute the polarization degree of the inverse Compton emission from 
isotropic electrons with arbitrary energy spectra scattering monoenergetic unidirectional 
photons. The authors note that the results imply a vanishing polarization degree for the 
inverse Compton emis s ion fr om unpolarized target photon beams. In a follow-up paper, 



Bonometto fc Saggion 



(119731 ) (referred to as BS in the following) describe a study of the 
polarization of the synchrotron self-Compton (SSC) emission from electrons in a magnetic 
field inverse Compton scattering their synchr otron emission. The ca l culati on involves the 



numerical integration of the BCS equations. iNagirner &: Poutanenl (jl993[ ) confirm the 



analytical results of BCS and cast the results into the form of a scattering matrix which 
operates on the Stokes vector of the target beam. The authors e yaluate th e scatt ering 



matrix for power law and Maxwellian electron distributions, and 



Poutanenl ( 119941 ) uses 



these results to discuss some aspects of the polarization of SSC emission. 

The shortcoming of all the analytical calculations is that it is not clear at which Lorentz 
factors the results start to be valid, and how accurate the predicted polarization degrees 
are in any specific case. Various authors have thus d eveloped codes to numerically evaluate 
the polarization degrees. iBegelman fc Sikoral ( 119871 ) describe a formalism to numerically 



compute the polarization of inverse Compton scattered unpolarized target beams in the 
Thomson regime. They use this formalism to study the polarization of electron beams 
("AGN jets") with finite opening angles and with different internal structures (e. g. filled 
and hollow cones). Most of their discussion focuses on situations in which the observer is 
either not looking directly into a modestly relativistic electron beam, or that the intensity 
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of the electron beam varies strongly over angular scales of 7 l . 



Celotti k Matt 



199 



) (called CM in the following) extend the formalism of 
Begelman fc Sikoral (jl987l ) to treat partially polarized beams and present numerical 
results concerning the polarization of SSC emission from power law electron populations. 
They find that SSC emission is highly polarized - although somewhat less than the 
synchrotron emission. The polarization direction of the SSC emission aligns with that of 
the synchrotron emission, i. e. the electric field vectors of the synchrotron and SSC emission 
are perpendicular to the projection of the magnetic field lines in the PF onto the plane of 
the sky. While the results depend somewhat on the electron Lorentz factors in the mildly 
relativistic regime, they are largely independent of 7 for 7>10. Their calculations also show 
that unpolarized synchrotron photons give rise to SSC emission with vanishing or very 
small (< 5%) polarization degrees for all but very low Lorentz factors. Although the results 
of CM deviate by ~ 15% from those of BS, the authors do not investigate the origin of this 



discrepancy. 



More recently, 



McNamara et al. 



McNamara et al. 



2008 



for additional 



((20091) (see 

information on the numerical algorithm) report on simulations of the emission from 
AGNs for specific geometries and for various unpolarized (accretion disk emission, cosmic 
microwave background (CMB)) and polarized (SSC emission) target photon fields. Among 
other cases, they compute the polarization of inverse Compton scattered unpolarized CMB 
emission. They report a polarization of ~20% in contradiction to the predictions of BCS of 
a vanishingly small net-polarization. 

In this paper, we re-evaluate the polarization properties of inverse Compton and SSC 
emission based on both analytical and numerical calculations. We describe for the first time 
detailed comparisons of analytical and numerical results. The comparisons show in which 
regimes the analytical approximations can be used. Our discussion focuses on establishing 
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the general properties of inverse Compton scattered emission. In this spirit, we discuss SSC 
emission as a special case of inverse Compton emission. Although we show some results 
for mildly relativistic electrons (7- values of 2 and 5), we focus the discussion on the regime 
7>10, where the polarization degree and direction have simple dependencies on the relative 
orientation of the electron momentum, target photon momentum, and target photon 
polarization vectors. We use the numerical simulations to investigate the discrepancies 
between the analytical and semi -analytical results of BC S and BS on the one hand, and 



the numerical results of CM and 



McNamara et al. 



(120091 ) on the other hand. Furthermore, 



we present for the first time calculations of the polarization degree of inverse Compton 
radiation emitted by isotropic electrons in the Klein-Nishina regime. 

In Section [2] we describe a generally applicable formalism that can be used to simulate 
inverse Compton emission from arbitrary electron and target photon distributions. The 
formalism can be used for arbitrary electron Lorentz factors and for scatterings in the 
Thomson and in the Klein-Nishina regimes, and permits the simulation of SSC processes. In 
Sections 13.11 and 13.21 we review the results from BCS and BS for the Thomson regime, and 
use them to derive analytical expressions for the intensity and polarization of the inverse 
Compton emission from monoenergetic electrons incident on monoenergetic unidirectional 
target photons. In Section 13.31 we perform detailed comparisons of the analytical and 
numerical results. The discussion covers the emission from monoenergetic electron and 
photon beams, from isotropic power law electrons scattering monoenergetic unidirectional 
photons, and from isotropic power law electrons scattering power law synchrotron emission. 
We scrutinize the polarization properties of inverse Compton scattered unpolarized CMB 
photons in Section I3.4[ and study the polarization properties of inverse Compton radiation 
emitted in the Klein-Nishina regime in Section 13.51 Section H] ends the paper with a 
summary and a discussion of the scientific potential of AGN observations with GEMS. 
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2. Monte Carlo simulations of inverse Compton processes in the Thomson 

and Klein-Nishina regimes 

2.1. General approach, Stokes vectors, and useful functions 

The numerical calculation uses a Monte Carlo approach that avoids some of the 
problems of earlier calculations based on integration over the entire phase space. Celotti & 
Matt (1994) use the latter approach and describe its difficulties: "For high values of the 
Lorentz factor 7, the functions in the integrals are strongly peaked in the electron direction. 
This requires the adoption of a carefully sampled angular grid over which to perform the 
angular grid...". In contrast, the Monte Carlo approach generates events with a relative 
frequency matching their natural occurrence and thus makes good use of the computational 
resources, even in the extremely relativistic regime. The simulation of inverse Compton 



processes resembles in some aspects the simulation of Compton scatterings (e.g. 



Matt et al. 



1996 



Depaola 



20031 ) with the added complication of a Dopper boost from the PF to the 
electron rest frame (EF) and back. The calculation begins and ends in the PF. We use a 
right-handed set of unit vectors C = {x, y, z} to define a reference coordinate system in the 
PF. Initially, we consider monoenergetic unidirectional electrons with an energy-momentum 
four- vector p = ( / ~fm e c, p) scattering a monoenergetic unidirectional beam of photons with 
the four- vector k = (co/c, k) with k = z (see Figured]). 

We use unit-less Stokes vectors s = (i, q, u) to track information about the statistical 
weight (i) and the linear polarization properties (via q/i and u/i) of the simulated photons. 
Conceptuall y, we adopt the q uantum mechanical treatment of the Stokes parameters 



described by lMcMasterl (jl96ll ). Accordingly, the Stokes parameters are expectation values of 
measurement results on statistical ensembles of incoherent photons. The additive property 
of the Stokes parameters makes them the tool of choice for the numerical treatment, as the 
results from many scattering processes can be combined by adding the Stokes vectors. 
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In the following, we will define Stokes vectors s a relative to various coordinate systems 
C a = {x a , y a , z a } with the z-axes always aligned with the photon propagation direction. 
When defining new coordinate systems below, we will only give the y and z unit vectors, as 
x a can always be computed from x a = y a x z a . We use the convention that a Stokes vector 
s a = (1, 1, 0) refers to a 100% linearly polarized photon beam with the electric field vector 
e parallel to q+ = y a . The Stokes vector s a = (1, —1, 0) refers to linearly polarization along 
q = x a (see Figure [2]). Correspondingly, s a = (1,0,1) and s a = (1,0,-1) refer to the 



electric field vector parall el to u + = 



Following the notation of 



(y a . ~ x a )/\/2 and u = (y a + x a )/\/2, respectively 



McMasterl (jl96ll ). the rotation of the Stokes vector s = (1, 1, 0) 



by an angle x clockwise looking along — k gives 

/l\ / 

si = M[ X ] 



1 

\0J 



1 \ 

cos 2x 
\ -sin2x ) 



(1) 



M[ X ] 



(2) 



where we introduced the matrix 

I 1 \ 

cos 2x sin 2% 
\ — sin 2x cos 2% / 
Equation (JTJ) shows that the angle x between the q + -direction and the polarization vector 
is given by the expression 



X 



arctan ■ 



— u 



Q 



(3) 



where some care has to be taken to select the right branches of the arctan-function to avoid 
jumps of the inferred x~ values. The equation for the polarization degree II reads 



n 



(4) 



From the above definitions, it follows that a Stokes vector s transforms as 



si = M[x] s 



(5) 
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for a coordinate transformation from Co to C\ resulting from a counterclockwise rotation 
around the z-axis of C (looking along -z ). 

For convenience, we introduce three functions. The first function operates on a vector and 
gives the normalized version of its argument 

A[a] = "4=1 (6) 
ya • a 

the second gives the component of one vector perpendicular to a second vector normalized 
to unit length 

P[o,b] = n[a- (a - b)b]; (7) 

the third function gives the angle between the y-axes of two coordinate systems which share 
a common z-axis: 

A[C a , C h ) = -sign(x a ■ y b ) arccos (y a ■ y b ). (8) 
A positive value indicates that Cb results from a counterclockwise rotation around z a . 



2.2. Simulation of one inverse Compton process 

For each photon, an inverse Compton scattering process is simulated with the following 
procedure. See Figure [3] for sketches of the scattering process in the plasma and electron 
frames. We dash all quantities in the EF. 

(1) A random number generator is used to draw the direction p (6 e ,<j) e ) of the scattering 
electron in the PF. We track the likelihood of an electron-photon interaction by multiplying 
the Stokes vector of the photon by the factor (l-/3/x kp ). Here, (3 is the velocity of the 
electron in units of the speed of light and /ik P = k • p is the cosine of the angle between the 
propagation directions of the target photon and the electron. 

(2) In the second step, we Lorentz transform the photon's wave vector and Stokes vector 
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into the EF using 

1 - Hkp/3 

and 



= (9) 



k' = /4 P 'P+ vi-Kv 2p [ k 'P]- ( 10 ) 

The photon frequency transforms as 

u/ = 7(1 - /3/ikpV (H) 

with 7 = ^1/(1 -/3 2 ) being the Lorentz factor of the electron. 

De Young showed that Lorentz boosts do not change the linear and circular polarization 
properti es of a statistica 



degrees ( iDe Young 



ensemble of photons including the linear and circular polarization 
19661 ). As we use a Stokes vector to track the statistical weight of a 
scattering process and information about the polarization direction, but not to track a 
beam intensity (radiation power per solid angle, area, and time), the transformed Stokes 
parameter i is not affected by the Lorentz boosts. We transform the Stokes vector into the 
coordinate system C*k p = {n[k x p],k} with the y-axis perpendicular to the (k-p)-plane: 

s kp = M[0 C ] s. (12) 

The Stokes vector s^p in the PF is identical to the Stokes vector s^'p defined in the EF with 
regards to the coordinate system Ck' p = {n[k x p], k'}: 

S k'p = S kp- (13) 

This important result follows from the Lorentz invariance of the polarization degree and 
from symmetry considerations alone. The latter imply that a beam polarized in the 
(k-p)-plane in the PF will be polarized in the same plane in the EF, because there is no 
preferred polarization direction perpendicular to the plane. Thus, the invariance of the 
polarization degree implies that a beam with q = —x and u = in the PF will have q = —x 
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and u = in the EF, i.e. that q is Lorentz invariant. Th e Lorentz invariance of II and q 



19871 ) imply 



implies the Lorentz invariance of u. Equations (6a-c) of (IBegelman fa Sikoral 
the same result, but were given without justification. 

(3) The third step implements the Compton scattering process in the EF. We denote the 
wave vector of the scattered photon with V, and use a random number generator to draw a 
random direction for the spatial part of I' assuming a constant probability per solid angle. 
With fj,y V = k' • 1', the scattering angle is 9' s = arccos [i' k , v . We use Compton's scattering 
formula A A = ^^(1 — ^° calculate the frequency of the scattered photon: 

1 



e' 1 



AVi' 



-co 



(14) 



where e' = hu'/m e c 2 is the target photon energy in the EF in units of the electron rest 



mass. We can now use the part of Fano's Compton scattering mat rix relevant for 



polar i zation to calculate the Stokes vector of the scattered photon (jMcMasterlll961 



mear 



Fano 



1949 



1957J): 



s[ v = F[9' s ,e',e' s ] M[A[C Vp ,C lv }} s' k , p . (15) 

The matrix M transforms the Stokes vector s^, into the coordinate frame CW = 
{n[l' x k],l'} with the y-axis perpendicular to the (k'-l')-plane - a prerequisite for using 
Fano's matrix. Fano's matrix F describes Thomson scattering in the Thomson and 
Klein-Nishina regimes and reads: 

/ l + cos 2 ^ + (e / -e^(l-cos^) 



\ 



sin 2 ^ 



V 



sin 2 ^ 
1 + cos 2 ^ 

2cos^ J 



(16) 



with 



hw' s /m e c 2 being the energy of the scattered photon in the EF in units of the 



electron rest mass. We did not include here the numerical factor r 2 /2 (with r being the 
classical electron radius) in the definition of F, as we are not interested in absolute fluxes in 
the following. Note that Fano's matrix accounts for the statistical weight of the scattering 
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process by modifying the i-component of the Stokes vector; the resulting Stokes vector is 
defined relative to Cyy. 

(4) The final step consists in back-transforming the results into the PF. The cosine of 
the angle between the scattered photon direction and the electron direction fx' Vp , —V ■ p' 
transforms as follows into the PF: 

The scattered photon direction in the PF is 



i = w+a/HVpP'p]- (is) 

The PF frequency of the scattered photon is 

u s = 7 (l + /Kp'K- (19) 

After inferring the Stokes vector of the scattered photon into the coordinate system 
CVp = {n[i' x p],l'} with the y-axis perpendicular to the (l'-p)-plane: 

s[, p = M[A[CV k ,,CV P ]] s|, k „ (20) 
we Lorentz transform the Stokes vector into the PF relative to C\ v = {n[l x p], 1} 

sip = s{, p . (21) 

Finally, we compute the Stokes vector relative to the coordinate system = {P[k, 1],1} 
with the y-axis aligned with the projection of the target photon's propagation direction in 
the sky: 

s Ik = M[A[C Ip ,C Ik ]] s lp . (22) 
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2.3. Simulation of SSC emission 



The simulation of SSC emission from electron and photon powerlaw distributions 
proceeds along similar lines. The calculations assume a uniform magnetic field B oriented 
at an angle $b to the line of sight. A random photon propagation direction k is drawn 
assuming a constant probability per solid angle. Subsequently, we generate a random photon 



frequency between u\ and U2 from a power law distribution with dN^/dE^ oc E- 



>+i) 



Similarly, an electron Lorentz factor is drawn from a power law distribution with 
dN e /d^ oc 7~ p . The initial Stokes vector is set to 

/ 1 \ 

s k B = (1 - sin a+1 6» B 



(23) 



n s 
V / 

with 8b being the angle between k and the magnetic field B and /z kp = k • p as above. The 
two multiplicative factors account for the relative likelihoo d of a electron-photon interact ion 



and for the scaling of the synchrotron emissivity with 8b (jGinzburg fc Syrovatskii 



1963 1 



We assume that the emission is linearly polarized along the q + -direction and the value of 
the g-parameter gives the polarization degree of synchrotron emission with a power law 
index a: 



IT 



1 + a 



(24) 



a + 5/3 

As synchrotron emission is polarized perpendicular to the B-field, s k B is defined with 
regards to the coordinate system C k B = {n[k x B], k}. The rest of the procedure is now 
equivalent to the one described in Steps (2)-(4) above with the modification that the Stokes 
vector Sk p is now given by 

s kp = M[A[C k B,C kp ]] SkB . (25) 

At the very end of the calculation, we compute the Stokes vector relative to a coordinate 
system Cm = {P[B, 1], 1} with the y-axis along the B-field direction in the PF projected 
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onto the sky: 

s 1B = M[A[C lp ,Q B ]] si p . (26) 

3. Analytical and numerical results 

3.1. Analytical results derived by BCS and BS 

We continue to use the same naming conventions as in the previous section. The 
frequency of an inverse Compton scattered photon can be computed from Equations ( fTTT) . 
(|14p . and ([T91 . In the Thomson regime e' — > we get 

oo s =7 2 (1-/W(1 + /3^ V V. ( 27 ) 

This equation gives the following minimum and maximum scattered photon frequencies as 
function of 7, u>, and /^k p : 

UVin = ^—i^u (28) 



and 

W m ax = 27 2 (1 - /3/i kp )w, (29) 

respectively. We will see below, that w max sets a useful scale to characterize the frequency 
dependence of the intensity and polarization degree of inverse Compton emission. 

From the last equation we infer that the minimum Lorentz factor 7 m ; n required to 
produce photons of frequency co s is for 7 3> 1 and j3 « 1: 



7 ™" = pzjr^-y (30) 

BCS and BS derived an expression for the intensity Jg s (energy emitted per unit time, 
unit volume and unit frequency interval) of a unidirectional monoenergetic photon beam 
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polarized along the vector e scattered by an isotropic electron population in the direction 1 
into a state of linear polarization defined by the polarization vector e s : 

(S! + S 2 ) + S 2 V (31) 



Je s K 7min 



* a , (k'6 s )(l-6) 



l-(k-l) 



with K = |/icr 2 cm 3 and 7 m i n as defined above. The factors Si and S 2 depend on the 
electron energy spectrum: 



Si = J m( 7 ) fx 2 - + 2^ rfx 



(32) 



S 2 = / 1 m( T ) (1 9 X)2 rfa:. (33) 



;r 2 



The function m depends on dn e /d , y, the differential electron spectrum per unit volume, as: 

m( 7 ) = ^ 7 " 2 , (34) 



and x is defined as 



x = ^. (35) 

7 



Summing Equation (!3~T|) over two orthogonal polarization directions gives the total intensity 
of the inverse Compton scattered radiation: 

Je s + Je s± = iT7min-(Si+3S 2 ). (36) 

BCS evaluated the linear polarization fraction for three exemplary cases with different 
initial polarization directions relative to the scattering plane (e s parallel, at 45°, and 
perpendicular to the scattering plane). The main results of the calculations are: 

(i) The intensity of the inverse Compton radiation is independent of the polarization degree 
and direction of the target photon beam, and does not depend on the scattering direction. 

(ii) The linear polarization degree of the inverse Compton radiation depends on the degree 
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but not the direction of the polarization of the target photons, and does not depend on the 
scattering direction. For a 100% polarized target beam, the linear polarization degree of 
the inverse Compton scattered radiation is 

"*» = (37) 

Note that although Si < is possible, the combinations Si + S2 and Si + 3S2 are equal or 
larger than zero. BCS evaluated IIbcs analytically for powerlaw electron distributions and 
numerically for the case of SSC emission. 

(hi) In contrast to Thomson scattering by stationary electrons, Inverse Compton scattering 
by non-thermal electron plasmas does not create polarization, and an initially unpolarized 
beam stays unpolarized after scattering. This result implies that if inverse Compton 
scattering reduces the polarization degree of a 100% polarized target beam to III, it will 
reduce the polarization degree of a target beam with a polarization degree of Il 2 to a 
net-polarization of 111 II2. 



3.2. Additional analytical results 

Equation (13TT) allows us to derive a closed expression for the polarization direction of 
the inverse Compton scattered radiation. The latter is given by the unit vector e s which 
maximizes the term in square brackets of Equation (13T|) with the constraining condition 
that § s _L 1. If we pull out the vector e s from the term in Equation (131 j) in square brackets, 
we see that the condition is equivalent to minimizing the angular distance between e s and 
the vector 

v = g + — -A^k. (38) 
l-(k-l) 
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The angular distance is minimized for the projection of v onto the plane perpendicular to 
1. The polarization vector is thus given by the expression 

e s = P[v,l]. (39) 

Using this expression, it can be shown in a somewhat tedious calculation that the 
polarization angle x = arccos (e B • P[k, 1]) equals the angle between the initial polarization 
direction and the projection of 1 onto the plane perpendicular to k: 

X = arccos(g • P[l,k]). (40) 

While Equation ( I3"9~j) gives the polarization direction, Equation ([411]) gives only the absolute 
value of x Du t not the sign of x- With regards to Figure [U the polarization angle x is 
7r/2 — (f) c . Note that the polarization direction does not depend on the observed frequency. 

The BCS equations allow us to derive the frequency dependence of the intensity and 
polarization degree of the inverse Compton emission from monoenergetic electrons. The 
differential energy spectrum of a monoenergetic electron beam normalized to one electron 



per cm 3 is 



^ = 5(7 - 7b) cm" 3 . (41) 



The integrands in Equations (1321) and (I3"3l only contribute for 



We get after some algebra 



and 



7min /An\ 

X = . (42) 

7o 



2 ^" ~\~ Xn 

El = XjL ~ 3 (43) 

7o 3 



E 2 = (izfOI. (u) 
£o7cT 
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Note that both parameters Si and E2 depend through xq (and thus through 7 m i n ) on the 
observed frequency u s . We obtain the total intensity 

4-\ 



and the polarization degree 



J = K u. 20-23^+2^ 
u 70 2 



, ^b-i . (46) 



1 - 2x 2 + 2x 

We now introduce the dimensionless variable y as the observed frequency in units of the 
maximum frequency w max (70) A*kp) which electrons of Lorentz factor 70 emit in a certain 
direction: 



The intensity scales with y as 



y = (47) 



1 - 2y + 2y 2 

■J oc y 5 • ( 48 ) 

7(T 



We infer a differential photon number per unit frequency interval of 

dN y dN 7 dy 3(1 - 2y + 2y 2 



P^/ = -j^-f- = V 9 ^ (49) 



'max 



where we introduced the normalization factor 3/(2 w max ) so that J^ max p 10 du s = 1 for 



The polarization degree as function of y is 

V 2 

n = -. (50) 

70 l-2y + 2y 2 K ' 

Equations ( 1491) and (1501 give the intensity and polarization fraction of the inverse Compton 
emission from monoenergetic electrons and can be used as integration kernel to predict 
the spectral shape and polarization properties of the inverse Compton emission from more 
complicated electron distributions. 
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3.3. Comparison of numerical and analytical results in the Thomson Regime 

It is instructive to look at the numerical results for a monoenergetic unidirectional 
photon beam scattering off a monoenergetic unidirectional electron beam. Figure H] shows 
the Stokes parameters i, q, and u and the polarization direction \ for an electron beam with 
6 C = 85°, e = 0° and 7 = 100 impinging on a photon beam with frequency u = 10 12 Hz 
incident along the z-axis (k = z), 100% polarized along the y-axis with an initial Stokes 
vector s = (1,1,0). As explained in the previous section, the polarization direction is 
measured clockwise from the direction of k in the sky when looking into the inverse 
Compton scattered beam. 

In the electron rest frame, the intensity of the Thomson scattered photons resembles 
the torus-shaped intensity emitted by a Hertzian dipole. After back-transformation into the 
PF, the intensity of the inverse Compton emission shows two poles with zero intensity at 
an angular distance of 2/7 from each other bracketing the high- intensity radiation emitted 
in the direction of the electron direction. The q and u distributions look complicated, but 
encode simple information: the inverse Compton emission from monoenergetic unidirectional 
electrons is everywhere polarized to 100%, and beyond the poles of the Hertzian dipole, 
the polarization direction exhibits a smooth rotation with a periodicity of 7r/2 around the 
maximum-intensity direction. The frequencies of the scattered photons (not shown here) 
are approximately evenly distributed between co min and w max (dN^/dojs ~ const). 

For an isotropic electron configuration, electrons with initial directions within a few 
7 _1 from the observed direction contribute i, q, and u-patterns as shown in Figure HJ It 
is clear that summing up positive i-values and positive and negative q and u-values will 
result in an overall polarization degree smaller than unity. Figure [5] shows the intensity, 
polarization degree and polarization direction for an isotropic electron distribution with all 
other parameters being the same as in the previous example. The angular distribution of 
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the inverse Compton emission (not shown here) is constant up to a dependence oc (1 — 0fiid) 
(with = k • 1) owing to the relative motion of the electrons and photons. Figure [5] 
shows the polarization degree and polarization direction for all scattering directions. The 
simulations show that the polarization degree is II = 0.5 independent of the scattering 
direction and that the polarization direction depends on k, e, and 1 as predicted by 
Equation (13"9|) . 

For mildly relativistic electrons, the polarization properties do depend on the scattering 
direction. As an example, Figure E] shows the same as the previous figure, but for 7 = 2. 
At an angle 9\ ~ 60° to the target photon beam, the polarization degree dips and the 
dependence of x 011 the azimuthal angle of 1 reverses. The main difference between the 
mildly relativistic and highly relativistic cases is that in the latter p and k are almost 
always nearly antiparallel while in the former, the photons approach the electrons from a 
wide range of directions. In the mildly relativistic case, the polarization direction of the 
incident photon is thus not "rotated" around an axis perpendicular to the p-1 plane. 

In the next step, we compare the frequency dependence of the intensity and polarization 
degree and direction with the theoretical predictions from Equations (1491) . (1301 . and (!38H39p . 
respectively. We consider the radiation scattered into the direction 1 with the polar 
coordinates (85°, 0°) for 7- values of 2, 5, 10, 20, and 100. Figure [7] shows that Equation (H9l) 
describes the frequency dependence extremely well except for mildly relativistic electrons 
with 7 = 2. Furthermore, Equation (150]) gives a good description of the polarization degree 
- again, except for 7 = 2. Additional simulations at higher 7-values (7 = 2,500, 12,500, and 
62,500) show intensity and polarization degree distributions identical to the one shown for 
7 =100. We also verified that the polarization direction (not shown here) is independent of 
frequency. 

The case of a powerlaw distribution of electrons (dN e /dj oc 7 _p from 71 to 72) 
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scattering off monoenergetic unidirectional photons is discussed in the paper of BCS. In 
the region where the emitted energy spectrum is a power-law, the polarization degree is 
predicted to be 

n , w (i+p)(3 + p) (51) 

Equations (l4"9j) and (jl9|) give a frequency dependent polarization degree of 



with 7x' = Max(7i, 7mi n ). The equation is valid for 7 min < 7 2 ; no emission is found at 
frequencies with 7 min > 72 . For 71 < 7 min and 72 — » 00, we recover Equation ( 151]) . 

Figure [S] compares the prediction from Equation (1521) with simulated results for 
71 = 10, 72 = 10, 000, and electron indices p between 1.01 and 4. The energy spectra of the 
inverse Compton emission indicate three regimes: a low- frequency regime with 7 m ; n < 71, 
a power law regime starting at 7^ >7i, and a high-frequency regime with 7 min — > j 2 - In 
the low-frequency regime the polarization degree is low as the emission is produced at 
low ^/-values. In the power-law regime, the polarization is independent of frequency and 
agrees well with Equations (15 ip and ( I52p . The softer the energy spectrum, the higher is 
the polarization degree. In the high-frequency regime, the polarization increases towards 
II = 1 as u s — > w max (72) because the emission is produced at y-values close to 1. The 
results validate the predictive power of Equation ( I5TI) for the power law regime, and of 
Equation (IS"!?]) for all frequencies. 

We now turn our attention to the polarization of SSC emission from an isotropic 
electron population emerged in an uniform magnetic field B. We assume that the 
electron and synchrotron number densities are well described by power law distributions 
dNe/d'j oc 7~ p and dN^/dcu oc oj~( a+1 ^ over the relevant energy and frequency ranges. 
The main objective of our discussion is to explain the polarization degree of the SSC 
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emission from the polarization properties of the synchrotron emission and the general 
results concerning the depolarization effect of inverse Compton processes described above. 

As discussed in Section I2.3[ the intensity of the synchrotron emission is proportional 
to sin a+1 #b, and its polarization degree is lis. In the power-law regime of the SSC 
radiation, the inverse Compton scattering reduces the polarization degree of synchrotron 
photons of a certain fixed frequency and propagation direction by the factor n p i. As 
the polarization direction of the inverse Compton emission depends on the propagation 
direction and polarization direction of the synchrotron photons, the net-polarization degree 
of the emerging SSC emission will be lower than ripills and we write 

n SS c~^n pl n s . (53) 

The factor rj depends on the angle $b between the magnetic field and the line of sight, 
because the magnetic field direction determines the intensity and polarization pattern of the 
synchrotron emission. Furthermore, i] depends on the spectral index a of the synchrotron 
emission because of Doppler boosting effects. 

We can derive an estimate for rj in the following way. Relative to the coordinate 
system C\b with the y-axis parallel to the B-field projected onto the plane of the sky and 
the z-axis pointing towards the observer, the w-parameter of the SSC radiation vanishes 
owing to the symmetry of the problem. The factor 77 is then proportional to the absolute 
value of the g-parameter and can be computed by integrating the g-parameter over the 
angular distribution of the synchrotron emission. Note that p is approximately parallel to 1 
for most of the contributing electrons. The expression for 77 reads: 

| / c% (1 - /^ kp )" +1 (sm ObY +1 cos 2 x I 

v /dn fc (i-/?// kp )«+ 1 (sm0 B ) Q+1 ■ 1 j 

The integrals run over the directions k of the synchrotron photons. The first two factors of 
the integrands are weighting factors. One power of (1 — /?/%,) accounts for the dependence 
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of the Compton scattering probability on the relative velocities of the photons and electrons. 
The additional power a stems from the fact that the synchrotron photons have an EF 
frequency ui' oc (1 — /3/ik P ) according to Equation (ITT)) , and the number of target photons 
available at a certain u' in the EF scales proportional to dN J^ ^ = dN<yL ^ " oc (1 — /3 /%,)"• 
The factor (sin^B) " 1 " 1 accounts for the dependence of the synchrotron emissivity on the 
angle 9-q between k and B. x is the polarization direction of the inverse Compton scattered 
emission computed with the help of Equation ( 1381) . and we made use of the identity 
q = cos2x from Equation ([1]). The denominator normalizes the expression. The numerical 
evaluation of the integral in the nominator gives negative values, corresponding to parallel 
synchrotron and SSC polarization directions - perpendicular to the projection of the 
magnetic field in the PF onto the sky. 

A frequency dependent estimate of the polarization degree can be derived from 
Equations flUD and (150)) : 

_ \q\ _ |/^/;;^r 7 ?rf 7o ^5(i-/3 / , kp )(sin^) Q+1 P7on s n 70 cos2 X 
UssciUs) " T " / dn t C du m d l0 (1 - /W (sin e B )^ P70 

(55) 

The integrals run over the directions of the synchrotron photons, the synchrotron photon 
frequencies and the electron Lorentz factors. Only electrons with Lorentz factors 7 > 7^ = 
Max(7!,7 min ) contribute and we set the value of the inner integral to zero if 7J > 72 . The 
first three factors of the integrand weight according to the electron and photon number 
densities and the interaction probability, the fourth factor is proportional to the synchrotron 
emissivity, and the fifth factor gives the photon density of the inverse Compton emission. 
In the integral of the nominator we multiply all these weighting factors with the fractional 
polarization of the synchrotron photons and the inverse Compton emission, and multiply 
with the contribution of the considered emission to the g-parameter at frequency u s . 

The numerical evaluation of the SSC polarization follows the description in Section [2731 
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Figure |9] presents the intensity and polarization degrees of the SSC emission for u\ = 10 9 Hz, 
u 2 = 10 13 Hz, 7! = 10, 72 = 10 5 , for (a = 0.5, p — 2) and (a = 1, p = 3), and various 
t?B-values. After low values at the lowest frequencies, the polarization degree of the SSC 
emission Hssc reaches a rather stable plateau with a gradual increase before an eventual 
peak at the highest frequencies. The polarization degrees calculated with Equation ( l55l) 
(Figure [HJ right panel, dotted lines) agree well with the results from the Monte Carlo 
simulations. The right panel also shows the SSC polarization degrees presented in BS and 
CM. Whereas the polarization degrees of BS agrees well with our calculations, those of CM 
deviate by up to ~15%. We explain the disagreement with the fact that CM measure the 
polarization degrees at frequencies where electrons with 7 < 10 still play an important role. 

Figure [10] compares the ^-dependence of the emitted radiation with the predictions 
of Equations (I53H55I) . The Monte Carlo simulations and Equation (153]) give frequency 
dependent polarization degrees, and we used here the polarization degrees averaged from 
10 16 to 10 18 Hz. The Mo nte Carlo simulations and Equations (15311551) give very consistent 
results. IPoutanenl (11 994 and private communication) suggested that r\ in Equation ( |53]) 
is given by 77 = sin 2 ^B- Our calculations show that the polarization degrees indeed scale 
proportional to sin 2 $b, however even for $b = vr/2 we find that 77 <1. The parameterization 
rissc = ^0 sm2 n pl ng with rjo = 0.79 for (a = 0.5, p = 2) and rj = 0.85 for (a = 1, p — 3) 
gives a good description of the results shown in Figure [101 

3.4. The polarization degree of the inverse Compton emission from 

unpolarized target photons 

One important implication of the analytical calculations of BCS and the numerical 
results shown in Figure [5] is that the polarization degree of inverse Compton emission off 
unpolarized radiation fields vanishes as long as the electron Lorentz factors are >10. The 



-26 - 



reason is that the scattering of two beams polarized perpendicular to each other will give 
two inverse Compton beams polarized perpendicular to each other, adding up to a beam 
with a vanishing net- polarization. We t ested this prediction explicitly for a case similar to 



the one discussed by 



McNamara et al. 



(120091 ) . We consider a jet of an AGN at redshift 



z — 2, scattering unpolarized CMB photons into the X-ray band. The jet bulk Lorentz 
factor is Tj = 5 and the jet is oriented at an angle of i?j = 80° towards the line of sight. The 
observed inverse Compton radiation leaves the jet at angle of 80° in the AGN frame (AF), 
and an angle of 166° in the PF. 

We generate an isotropic distribution of CMB photons in the AF with a frequency of 
(1 + z) 1.6 x 10 11 Hz = 4.8 x 10 11 Hz, and Lorentz transform the photon directions and 
photon frequencies into the PF. As the polarization degree is Lorentz invariant, the target 
photons are unpolarized in the PF, and we draw a random polarization direction in the 
PF. Subsequently, we simulate the Compton scattering taking into account the interaction 
probability as function of the electron and photon velocity vectors, and accumulate the 
Stokes parameters of the radiation emitted into the direction of the observer. In the last 
step, we Lorentz transform the frequency of the emerging radiation into the AF, and 
redshift it into the observer frame. 

We simulated 10 million events with an electron Lorentz factor of 7 = 10, 000. The 
choice of 7 produces inverse Compton emission with F v and z/i^-peaks in the 1-10 keV 
energy range. We calculate a net-polarization of the inverse Compton emission of 0.26%. 
Note that the polarization degree is positive definitive, so that even an arbitrary large 
number of simulations will always produce positive polarization degrees. We divided the 
simulated data set in 100 sub-sets, and verified that the q and u parameters of the sub-sets 
exhibited mean values statistically consistent with 0, as expected for zero net-polarization. 
These results do not confirm polarization degrees of ~20% of inverse Compton scattered 
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unpolarized photons as reported by 



McNamara et al 



(120091 ) . 



3.5. Inverse Compton emission in the Klein-Nishina regime 

The Klein-Nishina regime starts when the EF energy of the target photon is of the 
order or exceeds the electron rest mass energy: e'>l. In the deep Klein-Nishina regime 
(e' ^> 1), Equations (JIT]) . ( JT4j) . and ( JT9]) imply that head to head collisions of electrons 
and photons are likely to produce photons with an energy close to that of the incoming 
electrons. If we use the Fano Matrix with a Stokes vector of (1,1,0), we see that the 
intensity of the emission scattered into the direction of the scattering electron scales with e' 
as (1 + e' — e' s ) ~ e /_1 , and its fractional polarization scales as (1 + e' — e^) -1 « e' _1 . 

The third term of the (1, l)-element of the Fano matrix introduces unpolarized emission 
and an emission pattern which depends in the EF mainly on the polar scattering angle and 
not on the azimuthal scattering angle. 

Figure [11] shows maps of the intensity distribution and the polarization degree of a 
monoenergetic unidirectional photon beam with initial Stokes vector (1,1,0) scattered by 
monoenergetic unidirectional electrons. We chose e = 1/100 (~ 5.1 keV), k = (0,0,1), 
7 = 2,500 and p = (cos 85°, 0, sin 85°) which gives e' = 22.8. The intensity distribution is 
centrally peaked and decreases rather monotonically with the distance from the peak. The 
Hertzian dipole pattern is strongly suppressed. As discussed above, the polarization degree 
at the peak of the intensity distribution is approximately 1/e'. The polarization degree 
shows two dips with vanishing polarization degrees. These dips correspond to scattering 
directions along the polarization direction of the target photon in the EF. 

Figure [12] presents the intensity and polarization degree for a monoenergetic and 
unidirectional photon beam with an initial Stokes vector (1,1,0) scattering off monoenergetic 
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isotropic electrons. We choose e = 1/400 (~ 1.3 keV), k = (0,0,1), and an observer at 
1 = (cos 85°, 0, sin 85°). We simulate 7-values between 10 and 62,500 leading to e'-values 
(computed with /i^ as proxi for /x^p) between 0.02 and 143. All distributions are shown as 
function of the frequency u s in units of the maximum frequency allowed kinematically 

Wmax,KN — : , . • {00) 

1 + 47W 

Deeper into the Klein-Nishina regime (at larger e'-values), the intensity distribution is more 
and more peaked towards w mffi) KN, and the polarization is strongly suppressed. Figure 
[TBI displays the net-polarization degree as a function of 7. To guide the eye, the figure 
compares the numerical polarization degree results with the function II = 0.5/(1 + e'). The 
latter function agrees with the numerical results for scattering in the Thomson regime, but 
deviates considerably in the mild and deep Klein-Nishina regimes. 



4. Discussion: Application to GEMS observations 

In this paper, we describe a general formalism that can be used to study the polarization 
of inverse Compton emission, including SSC emission. The comparison of the numerical 
results with the analytical results of BCS and BS validate the analytical results for 7 > 10. 
The numerical simulations show that Equations (l4"5]L ( l50l) . ( 1551) . and ( J38ll3~9l can be used 
to predict the polarization properties of inverse Compton and SSC emission. An important 
qualitative result is that inverse Compton scattering of unpolarized target photons does not 
create a polarized signal as long as the electron distribution is isotropic over angular scales 
of 7 _1 and the electron Lorentz factors exceed minimum values of ~10. 

The GEMS X-ray polarimetry mission will have the opportunity to make major 
discoveri es concerning AGNs, includ ing constraints on the structure of accretion disk 



coronae (jSchnittman fc Krolik 



20101 ). the role and structure of the magnetic field in AGN 
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jets, and the origin of the low-energy an d high-energy components of the continuum 



emission from blazars (IKrawczynski 



20111 ). Most relevant in the context of this paper are 
GEMS observations of the spatially unresolved continuum emission from blazar jets, because 
blazars are very bright X-ray sources. Two classes of objects are of particular interest. Flat 
spectrum radio quasars (FSRQs) and low frequency and intermediate frequency peaked BL 
Lac objects (LBLs and IBLs) with energy spectra with a uF u -peak in the IR/optical/UV 



band and one in the MeV to GeV band, and high frequency peaked BL Lac ob 
with one ^i^-peak in the UV/X-ray band and one in the GeV/TeV band (see 



ects (HBLs 



Abdo et al 



2010 



Giommi et al 



2002 



for compilations of spectral energy distributions). 



In the case of HBLs (e.g. the X-ray bright sources Mrk 421, Mrk 501, 1ES1959+650, 
PKS 2155-314, and PKS 1218+304), GEMS will sample the low-energy emission 
component, presumably of synchrotron origin. Optical blazar observations exhibit rather 



high polarizat ion degrees, often c 



emission (e.g. 



Angel fc Stockman 



ose to the theoretical maxim um lis f° r synchrotron 



1980 



Scarpa &: Falomo 



19971 ) . If X-ray emission from 



HBLs is also synchrotron emission, a commonly made assumption, GEMS observations 
should reveal X-ray polarization degrees and directions similar to those observed in the 
optical band. Indeed, the polarization degrees could be even higher, as X-ray emitting 
electrons radiatively loose their energy on shorter time scales than optically emitting 
electrons. As the X-ray emitting electrons have thus less time to travel away from the 
acceleration sites, the X-ray bright regions are expected to be smaller than the optically 
bright regions. The more uniform magnetic field in the smaller X-ray bright regions should 
lead to a higher polarization degree of the synchrotron emission. The same reasoning 
suggests the possibility that continuous swings of the polarization direction as occasionally 



observed in the optical (e.g. 



Marscher et al. 



20081 ). might be observed more often in the 



X-ray band. The polarization swings have been interpreted to arise from the movement o f 



a helical magnetic field which threads the jet through stationary shocks ( IMarscher et al 
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20081 ). Thus, GEMS has the potential to deliver observational evidence for a helical 



magnetic field at the bases of jet s, a predictio n of magnetic models of jet formation, 



acceleration, and collimation (e.g. 



Spruit 



20101 . and references therein). 



For FSRQs (e.g. the bright sources 3C 279, S 52116+81, and 1ES 0836+710), LBLs 
(e.g. the bright sources BL Lac, ON 231, and OJ 287), and some IBLs, GEMS will sample 
the high-energy vF v - component, presumably of inverse Compton origin, and can be used 
to constrain the origin of this component. If the emission is of SSC origin, the inverse 
Compton emission should track the polarization degree and polarization direction of the 
synchrotron emission. Figures 151 and [TU1 show that the SSC emission can be highly polarized. 
In external- Compton (EC) models, the dominant target photons come from the accretion 
disk, from the corona, from the broad line region clouds (BLR), or the CMB. As discussed 



in Section 13.31 unpolarized target photons would result in very low (<C 1%) polarization 
degrees. Although reflection off the BLR clouds would polarize the target photons coming 
from a particular direction, averaging over an axisymmetric cloud configuration would 
lead to a low net-polarization of the target photons and thus to a rather low polarization 
degree of the EC emission. The detection of high polarization degrees would thus strongly 
favor a SSC origin of the X-ray emission. However, the detection of vanishing or low 
polarization degrees would not exclude the SSC scenario, as the SSC process can produce 
small polarization degrees if the magnetic field is aligned with the line of sight. The SSC 
hypothesis would even be consistent with low polarization degrees observed for many 
sources, as it cannot be excluded that the blazar-zone magnetic field is preferentially 
aligned with the line of sight. If the polarization degree is low but non-zero, simultaneous 
optical/X-ray measurements of the polarization degree and direction should be able to 
decide between an SSC and an EC origin, as the former predicts a good correlation while 
the latter does not. 
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In the discussion so far, we have assumed that the magnetic field, and the target 
photon and electron distributions do not change over the emission region. If information 
about the spatial (and/or temporal) distributions of these quantities is available, the 
polarization degree can be calculated with a calculation similar to Equation ( |55|) . In this 
case, one needs to integrate the expression in the nominator (inside the absolute-value 
brackets) and in the denominator over the emission volume. Inside the integrals, the 
magnetic field as well as the target photon and electron densities should be evaluated 
at the retarded times. As Equation f )55|) gives only \q\/i, one has to calculate also \u\/i 
by replacing cos 2% in the nominator by sin 2%. The polarization direction and degree 
can then be computed with Equations ([3]) and (jl]), respectively. In most spatial 
variation of the source properties will lead to a large number of free model parameters, 
and will render the unambiguous interpretation of observational data more difficult. 
Some of the most interesting GEMS results may come from observations of polarization 
degrees close to the maximum theoretically possible polarization degrees, which indicate a 
uniform emission region and make it possible to constrain models effectively. The reader 
interested in non-uniform emission regions should refer to the related discuss ions of the 
polar i zation of synchrotron emission from non-uniform emission regions (e.g. 



1993 



Ruszkowski fc Begelman 



2002 



Begelman 



and references therein). 



The Monte Carlo approach presented in Section [2] can also be used to explore the 
polarization properties of inverse Compton emission in the case that the electron phase space 
distribution is not isotropic in any frame of reference. Future research could combine the 



the structure of blazar jets (see e.g. 



McKinnev & Blandford 


2009; 


Meliani & Keppens 


2009 



and references therein) with the results 



shocks (see e.g. 



Sironi &: Spitkovsky 



2009 



rom p article-in-cell simulations of relativistic 



2011 



and references therein) to derive estimates 



of the electron phase space distribution. Monte Carlo simulations as described in this paper 
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can then be used to derive observational signatures which can be compared to experimental 
data. 

Acknowledgments: HK thanks the members of the GEMS science working group, 
especially M. Baring, J. Poutanen, and J. Scargle for discussions on the polarization degree 
of synchrotron and inverse Compton emission. The author acknowledges helpful comments 
of an anonymous referee, and thanks A. Parvulescu and J. W. Krawczynski for proofreading 
the manuscript. He acknowledges NASA for support from the APRA program under the 
grant NNX10AJ56G, the DOE for support from its high-energy physics division, and 
support from the McDonnell Center for the Space Sciences at Washington University. 



Facilities: GEMS 



-33 - 

REFERENCES 

Abdo, A. A., et al. 2010, ApJ, 716, 30 
Angel, J. R. P., Stockman, H. S. 1980, ARA&A, 8, 321 
Begelman, M. C. 1993, Lect. Not.in Phys., 421, 145 
Begelman, M. C, & Sikora, M. 1987, ApJ, 322, 650 

Bellazini, R., Costa, E., Matt, G., et al. 2010, "X-ray Polarimetry: A New Window in 
Astrophysics", Cambridge University Press. 

Bjornsson, C.-L, & Blumenthal, G. R. 1982, ApJ, 259, 805 

Bjornsson, C.-I. 1985, MNRAS, 216, 241 

Black, J. K., Deines- Jones, P., Hill, J. E., et al. 2010, "The GEMS photoelectric x-ray 
polarimeters" , Space Telescopes and Instrumentation 2010: Ultraviolet to Gamma 
Ray. Edited by Arnaud, Monique; Murray, Stephen S.; Takahashi, Tadayuki. 
Proceedings of the SPIE, 7732, 77320X 

Blandford, R. D., Payne, D. G. 1982, MNRAS, 199, 883 

Bonometto, S., Cazzola, P., & Saggion, A. (BCS) 1970, A&A, 7, 292 

Bonometto, S., & Saggion, A. (BS) 1973, A&A, 23, 9 

Celotti, A., & Matt, G. (CM) 1994, MNRAS, 268, 451 

De Young, D. S. 1966, Journal of Mathematical Physics, 7, 1916 

Dean, A. J., Clark, D. J., Stephen, J. B., et al. 2008, Science, 321, 1183 

Depaola, G. O. 2003, NIMA, 512, 619 



-34 - 



Dolan, J. F. 1967, Space Sci. Rev., 6, 579 
Fano, U. 1949, J. Opt. Soc. Am. 39, 859 
Fano, U. 1957, Revs. Modern Phys. 29, 74 

Forot, M., Laurent, P., Grenier, I. A., et al. 2008, ApJ, 688, L29 
Ginzburg, V. L., Syrovatskii, S. I. 1965, ARA&A, 3, 297 

Giommi,P.,Capalbi,M.,Fiocchi,M., Memola, E., Perri, M., Piranomonte, S., Rebecchi, 
S., Massaro, E. 2002, A Catalog of 157 X-ray Spectra and 84 Spectral Energy 
Distributions of Blazars observed with BeppoSAX, Procs.: Blazar Astrophysics 
with BeppoSAX and other Observatories (Frascati, December 2001), eds. Paolo 
Giommi, Enrico Massaro and Giorgio Palumbo [arXiv:astro-ph/0209596], http: 
/ /www. asdc.asi.it/blazars/ 

Kaaret, P., Swank, J., Jahoda, K., et al. 2009, "The Gravity and Extreme Magnetism Small 
Explorer (GEMS)", Publication: Chandra's First Decade of Discovery, Proceedings 
of the conference held 22-25 September, 2009 in Boston, MA. Edited by Scott Wolk, 
Antonella Fruscione, and Douglas Swartz, abstract #125 

Krawczynski, H. 2011, Astroparticle Physics, 34, 784 

Krawczynski, H., Garson, A., Guo, Q., Baring, M. G., Ghosh, P., Beilicke, M., & Lee, K. 
2011, Astroparticle Physics, 34, 550 

Laing, R. A., 1980, MNRAS, 193, 439 

Laurent, P., Rodriguez, J., Wilms, J., Cadolle Bel, M., Pottschmidt, K., Grinberg, V. 2011, 
Science Express, Science DOI: 10.1126/science.l200848 

Lei, F., Dean, A. J., Hills, G. L. 1997, Space Sci. Rev., 82, 309 



-35 - 



Lightman, A. P., Shapiro, S. L. 1975, ApJ, 198, L73 
Lightman, A. P., Shapiro, S. L. 1976, ApJ, 203, 701 
Marscher, A. P., et al. 2008, Nature, 452, 966 

Matt, G., Feroci, M., Rapisarda, M., & Costa, E. 1996, Radiation Physics and Chemistry, 
48, 403 

McKinney, J. C, & Blandford, R. D. 2009, MNRAS, 394, L126 
McMaster, W. H. 1961, Reviews of Modern Physics, 33, 8 
McNamara, A. L., Kuncic, Z., & Wu, K. 2008, MNRAS, 386, 2167 
McNamara, A. L., Kuncic, Z., & Wu, K. 2009, MNRAS, 395, 1507 
Meliani, Z. and Keppens, R. 2009, ApJ, 705, 1594 
Meszaros, P., Novick, R., Szentgyorgyi, A., et al. 1988, ApJ, 324, 1056 
Nagirner, D. I., k Poutanen, J. 1993, A&A, 275, 325 

Novick, R., Kestenbaum, H. L., Long, K. S., et al. 1977, "OSO-8 X-ray polarimeter and 
Bragg crystal spectrometer observations", In: Highlights of astronomy. Volume 4 
- International Astronomical Union, General Assembly, 16th, Grenoble, France, 
August 24-September 2, 1976, Proceedings. Part 1. (A78-16417 04-88) Dordrecht, D. 
Reidel Publishing Co., 1977, 95 

Poutanen, J. 1994, ApJS, 92, 607 

Rees, M. J. 1975, MNRAS, 171, 457 

Rybicki, G. B., Lightman, A. P. 1986, Radiative Processes in Astrophysics, Wiley- VCH 



-36 - 

Ruszkowski, M., & Begelman, M. C. 2002, ApJ, 573, 485 
Scarpa, R., Falomo, R. 1997, A&A, 325, 109 
Schnittman, J. D., Krolik, J. H. 2010, ApJ, 712, 908 
Sironi, L., & Spitkovsky, A. 2009, ApJ, 698, 1523 
Sironi, L., & Spitkovsky, A. 2011, ApJ, 726, 75 

Spruit, H. C. 2010, "Lecture Notes in Physics", Berlin Springer Verlag, 794, 233 
Vlahakis, N., Kdnigl, A. 2004, ApJ, 605, 656 

Weisskopf, M. C, Silver, E. H., Kestenbaum, H. L., et al. 1978, ApJ, 220, L117 

Weisskopf, M. C, Eisner, R. F., Hanna, D., et al. 2006, "The prospects for X-ray 

polarimetry and its potential use for understanding neutron stars" , Paper presented 



at the 363rd Heraeus Seminar in Bad Honnef, Germany, arXiv:astro-ph/0611483 



This manuscript was prepared with the AAS IATj^X macros v5.2. 



-37- 




Fig. 1. — Sketch of the simulation set-up in the plasma frame. An unidirectional beam 
of electrons propagating along the direction p with polar coordinates (6 e , <f> e ) scatters an 
unidirectional photon beam propagating along the k = z. 



Fig. 2. — We define the polarization direction of a photon propagating along k with reference 
to a right handed coordinate system C = {x, y,z} with z = k. The electric field vector of 
photons with a Stokes parameter q/i =1 is aligned with the y-axis, that of photons with 
q/i =-1 is aligned with x-axis. Looking along -k, the directions with u/i—1 and u/i=-l are 
rotated counterclockwise by tt/A relative to the q/i =1 and q/i =-1 directions, respectively. 
The polarization angle x is measured clockwise from the direction with q/i =1. 



Fig. 3. — The Compton scattering process as seen in the plasma frame (PF) and in the 
electron (EF) frame. The PF electron and photon directions before scattering are p and k, 
respectively. The EF photon directions of the target and scattered photons are k' and 1', 
respectively. The direction of the scattered photon in the PF is denoted with 1. 
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Fig. 4. — Distribution of the Stokes i, q and it parameters and the polarization direction \ of 
the inverse Compton emission from a monoenergetic unidirectional electron beam (7 = 100, 
(6* e , e ) = (85°, 0)) scattering a monoenergetic unidirectional target photon beam propagating 
along the z-axis and polarized along the y-axis (co = 10 12 Hz, k = (0, 0, 1), s = (1, 1, 0)) as 
function of the polar coordinates of the scattered photon direction 1. x is measured relative 
to q + aligned with the projection of k onto the sky. 
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Fig. 5. — Polarization degree II and polarization direction x of the inverse Compton emission 
from a monoenergetic isotropic electron distribution (7 = 100) scattering a monoenergetic 
unidirectional target photon beam (to = 10 12 Hz, k =(0,0,1), s = (1, 1, 0)) as function of the 
polar coordinates of the scattered photon direction 1. \ is measured relative to q + aligned 
with the projection of k onto the sky. 
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Fig. 6. — Same as Figure EJ but for electrons with Lorentz factor 7 = 2. 
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Fig. 7. — Frequency dependence of the intensity (left panel) and polarization degree (right 
panel) of the inverse Compton emission from a monoenergetic isotropic electron population 
scattering a monoenergetic unidirectional photon beam. All parameter values are the same 
as in Figure HI except 7 = 2 (dashed red lines), 7 = 5 (dot-dashed blue lines), 7 = 10 (solid 
green lines), 7 = 20 (solid magenta lines), and 7 = 100 (solid orange lines). The dotted black 
lines display the analytical results from Equations ( 1491) and (!50l) . The frequencies are given 
as y = w s /w max in units of the maximum frequency for Thomson processes. The results are 
shown for photons emitting along the direction (6\, = (85°, 0). The plotting range of the 
right panel is restricted to the range where the statistical errors on the polarization degree 
are sufficiently small (a color version of this figure is available in the online journal). 
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Fig. 8.— Intensity per logarithmic energy interval (left panel) and polarization degree 
(right panel) of the inverse Compton emission from isotropic power law electron populations 
{dN B /drj oc 7~ p from 71 = 10 to 72 = 10 4 ) scattering a monoenergetic unidirectional target 
photon beam (u = 10 12 Hz, k =(0,0,1), s=(l,l,0)) into the direction (0 h fa) = (85°, 0). The 
thick to thin solid lines display the simulation results for electron spectral indices of p = 
1.01, 1.5, 2, 3, 4. In the right panel, the dotted lines show the semi-analytical results from 
Equation (|52|) . 
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Fig. 9. — Intensity per logarithmic energy interval (left panel) and polarization degree (right 
panel) of the SSC emission from isotropic power law electron populations (dN c /d , ~f oc 7~ p 
from 71 = 10 to 72 = 10 5 ) scattering synchrotron emission (dN^/dco oc cu~ ( - a+1 ' 1 from uj\ = 
10 9 Hz to u>2 = 10 13 Hz) from a B-field at 85° to the line of sight. Numerical results are 
reported for a = 0.5 and p = 2 (circles), and a = 1 and p = 3 (squares). The dotted lines in 
the right panel present the semi-analytical results from Equation ( 155)) . The dashed lines are 
the semi-analytical results from BS and the dashed-dotted lines report the numerical results 
from CM. Both papers do not give frequency resolved results. 
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Fig. 10. — Polarization degree of the SSC emission as function of the angle $b between the 
magnetic field and the line of sight for a = 1 and p = 3 (squares) and a = 0.5 and p = 2 
(circles), evaluated between 10 16 and 10 18 Hz. The dashed line displays the expectation from 
Equation (1551) . and the dotted line shows the results from Equations (153"]) and (1541) . 
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Fig. 11. — Stokes % parameter (left panel) and polarization degree II (right panel) of the in- 
verse Compton emission from a monoenergetic unidirectional electron beam (7 = 2, 500, 
(6* e ,0e) = (85°, 0)) scattering a monoenergetic unidirectional photon beam (e = 0.01, 
k =(0,0,1), s = (1, 1,0)) in the Klein-Nishina regime (e' =22.8). 
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Fig. 12. — Intensity (left panel) and polarization degree (right panel) of the inverse Compton 
emission in the direction (6\,<f)i) = (85°, 0) from monoenergetic isotropic electrons scattering 
monoenergetic unidirectional photons (e = 0.0025, k =(0,0,1), s = (1,1,0)). From left to 
right, the lines display the results for 7-values of 10, 100, 500, 2,500, 12,500, and 62,500. 
In the EF, the target photon energy in units of the electron mass is 0.02, 0.2, 1.1, 5.7, 29, 
143 for the six simulated cases. In the right panel, the frequencies are given in units of 
the maximum kinematically allowed frequency y = w s /w ma x,KN) and the plotting ranges are 
restricted to the ranges where the statistical errors on the polarization degree are sufficiently 
small. 
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Fig. 13. — Net-polarization degree of all emitted photons as function of the electron Lorentz 
factor (solid line and circles) for the simulations of Figure [T21 with additional simulations for 
7- values of 3.125 x 10 5 and 1.5625 x 10 6 . The dashed line displays the function II = 0.5/(l+e') 
showing that the net-polarization degree goes down approximately as the inverse of the EF 
target photon energy e'. 



